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Abstract Entanglement and non-locality are non-classical global characteristics of quantum states important 
to the foundations of quantum mechanics. Recent investigations have shown that environmental noise, even 
when it is entirely local in influence, can destroy both of these properties in finite time despite giving rise to 
full quantum state decoherence only in the infinite time limit. These investigations, which have been carried 
out in a range of theoretical and experimental situations, are reviewed here. 
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1 Introduction 

Entanglement and non-locality, which long have been of interest in the foundations of quantum theory, have 
been of renewed interest during the last two decades as the field of quantum information science emerged and 
matured. The techniques most often used in the latter context have recently been applied in the investigation 
of these properties, often with an eye toward the foundations of quantum theory. The investigation of the 
effects of noise on entanglement and non-locahty is a potentially important example of this, because it could 
prove more significant for the foundations of quantum theory than the study of the simpler related effect of 
quantum decoherence, from which they have been shown during its course to have important differences. 

Entanglement and non-locaUty are distinct properties that are often conflated. Therefore, it is important 
to clarify their relationship as well as their relationships to quantum coherence. Quantum coherence may be 
destroyed abruptly, for example, by strong projective measurements or asymptotically in time, for example, 
due to weak noise influences as a quantum system interacts with its enviroimient. It has been shown that even 
in the latter case, in which coherence disappears only in the infinite-time limit, entanglement and non-locality 
may still abruptly and completely disappear in finite time, phenomena named Entanglement Sudden Death 
(ESD) and Bell non-locality Sudden Death (BNSD), respectively. 

After their relatively recent discovery, ESD and BNSD have been explored in a variety of contexts, the- 
oretically and experimentally, in continuous and discrete systems, and when systems are subject to noise, 
both classical and quantum in origin, which may affect both amplitude and phase. The physical systems that 
have been examined are diverse and include electrons on a solid state lattice subject to electromagnetic fields, 
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photons traveling among collections of beamsplitters and mirrors, and atoms confined in electromagnetic cav- 
ities. The results suggest that ESD and BNSD are ubiquitous phenomena which must be kept in mind in any 
realistic physical analysis. 

This article provides a comprehensive survey of the investigation of ESD and BNSD. It is structured as 
follows. In Sec.|2] we briefly outline the relevant details of the theory of quantum coherence and open quantum 
systems in order to properly contextualize ESD and BNSD. In Sec.|3] we survey work concerning ESD in both 
theoretical and experimental contexts. In Sec.|4] we turn attention to the investigation of BNSD and highlight 
how similar and distinct from ESD it is. In Sec.|5] we conclude with an examination of how ESD and BNSD 
may be of important to the interpretations and foundations of quantum mechanics and quantum information 
science. 

2 Coherence 

In recent decades, quantum coherence has most often been examined from the perspective of its destruction. 
Decoherence has been shown to result in the degradation of quantum properties, corresponding to a transi- 
tion from the quantum domain to the classical domain. In general, decoherence may be classified into two 
main categories. The study of the first, intrinsic decoherence, follows an approach that modifies the standard 
quantum dynamics, leading to self-induced decoherence (SID); a stochastic term is added to the Schrodinger 
equation to induce spontaneous state-function collapse. Most studies of SID examine the following dephasing 
master equation for a state p and Hamiltonian H 1144111001 171. 

T being a theory-dependent parameter. In the energy basis, the solution is given by 

P(0= I =p„ye-'(^"-^"')'/''exp 
«,«'=o 

with p„_„'=o = {E„\p{0)\E'^), 7= T. The state purity is accordingly 

^(f)=trp2= £ |p„,„fexp 

n,n'— 

Provided that the temporal parameter t is a constant, the degradation of state purity, and so state coherence, 
occurs exponentially in time. 

The study of extrinsic decoherence, by contrast, uses standard quantum dynamics and considers the en- 
vironment surrounding the system of interest in interaction with it; the environment is typically modeled as 
a system of much larger dimensionality than that of the system of interest. Numerous classical and quantum 
correlations arise between the system and environment by unitary evolution. As a result, the reduced system 
state obtained by tracing over the environmental state generally exhibits decoherence and may also exhibit 
ESD or BNSD. 

2.1 Open quantum systems and noise 

The models discussed in this review, with a few exceptions, are primarily based on open quantum systems 
theory, which has as a core assumption that quantum systems are not viewed in isolation, but only with 
respect to and within some environment, that is, following the extrinsic decoherence approach. Not only is 
the environment's interaction with the system inevitable, it is not confined to a perturbative influence; the 
environment may strongly affect the system's time evolution and determine which quantities are observable. 
Even if the exact evolution equations of the quantum system coupled to the environment were solved, their 
complexity would preclude detailed and meaningful analysis in all but the simplest systems. 

Generally, open quantum systems theory considers a quantum system of interest, Psys, together with its 
environment, penv, in a joint state ptot = psys ® Penv, that evolves unitarily and independently according to their 
own internal Hamiltonians //sys and //env, respectively, with an interaction Hamiltonian //int, of the form 
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(3) 



^tot = //sys ® I + 1 ® /^env + ^^int 



(4) 
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The statistical density matrix of the joint system pst(0 evolves unitarily, Pst(0 = U{t)piotU^{0), with the 
unitary operator U{t) given by the Schrodinger solution, U (t) = exp [—ijQdt'H{t')\. During the interval 
[0,t), classical and quantum correlations build up between the system and the environment, with its much 
larger number degrees of freedom, and spread throughout the joint Hilbert space. The time-evolved density 
operator averaged over noise fields is p (t) = (Pst(f))noise- When attention is confined to the system of interest 
by tracing out the environment after unitary evolution, Psys = trenvPtot, the influence of the environment is 
seen as noise that causes dephasing and amplitude damping, asymptotically leading to complete decoherence. 

For simplicity and without loss of generality, this process may be modeled more concretely for purposes 
of calculation for specific physical situations in terms of the operator sum decomposition. The time-evolved 
density matrix is given by the following completely positive and trace preserving (CPTP) map, 

p{t)=^ [p (0)] = f (f ) p (0) (f ) , (5) 
/j=i 

where the operators in the decomposition Efi satisfy the respective positivity and trace preserving relations, 

Y.^El^{t)Efi{t) = I and Y.pE ^{t)E'^^{t) = I, f90l'65|. Provided that the specific forms of the operators £^ 
satisfy the CPTP relations, a set of non-unique operators that model a variety of quantum noise effects, such 
as dephasing and amplitude damping, as well as more general quantum evolution, can be constructed. The 
exact forms of the sets of operators {i?^ (f ) } are provided here if needed for clarity. An informative overview 
of this open quantum systems approach has been given by Benatti and Floreanini 1 18 |. 

Open quantum systems theory has been central to a better understanding of decoherence. This approach is 
motivated by the fact that environmental effects cannot always be treated as perturbative corrections. Despite 
environmental noise causing a quantum system to decohere only asymptotically as t ^ this noise may in 
finite time destroy uniquely non-classical properties such as entanglement or non-locality, that is, produce 
BSD or BNSD respectively. 



2.2 Coherence measures 

There exist a number of useful indicators of quantum coherence, both quantitative measures such as state 
purity and fidelity as well as qualitative methods such as the mere existence of off-diagonal terms in the density 
matrix. Those that are frequently used in the study of decoherence and disentanglement are summarized in 
this section. 

The purity of a state p is an easily computable measure of quantum coherence, which is given by 

,2 



= trp^ (6) 

For a t/-dimensional system, 1/d < 3^{p) < 1 ; the lower bound 1/dis reached only for the completely mixed 
state exhibiting no coherence, whereas the upper bound 1 is saturated by any pure state. As a state decoheres, 
its purity is reduced. The greatest degree of of decoherence occurs when a pure state, with initial unit purity 1 
later reaches purity 1 /d. 

The fidelity coherence measure for arbitrary mixed state density matrices pi and p2, is given by 



F(Pl,P2) = 



tr \ VP2P1VP2 



1 2 



(7) 



with < F(pi , P2) < 1 . The upper bound 1 indicates that pi and P2 are indistinguishable and the lower bound 
indicates that pi and p2 are orthogonal; the fidelity may be used in certain circumstances to find the time of 
the loss of coherence of a state from to its initial state by taking pi = p{t = 0) and p2 = p{t) represent the 
initial state and the later state, respectively. 

Although these measures are often suitable for quantifying the amount of quantum state coherence present, 
in many cases the concern is only whether quantum coherence exists or not; for that purpose, it is sufficient 
to observe the off-diagonal elements, only when they are all zero is there no state coherence in the particular 
basis in which they are represented. 
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3 Entanglement sudden death 



Let us now consider work specifically focusing on Entanglement Sudden Death (ESD). Sec. 3.1 briefly re 



views required background, including the most often used entanglement measures and the classification of 



quantum states. Sec. 3.2 discusses the initial discoveries of ESD in both continuous and discrete variable 



systems showing that decoherence and disentanglement are distinct phenomena in that they proceed at qual- 
itatively different rates. Sec. |3.3| surveys the subsequent theoretical development of the study of ESD. Sec. 
|3.3.7| discusses ESD as found in a range of different physical systems. Finally, experimental confirmation of 
ESD and related experimental issues are described in|3.4| 



3.1 Preliminaries 

In order to investigate Entanglement Sudden Death (ESD), one must bring to bear appropriate measures of 
entanglement, which are described in this section after introducing the relevant entanglement clases. 

3.1.1 Entanglement classes 

Because entanglement is a global property of a state and non-increasing under local operations, there are 
natural state classifications provided by state behavior under local operations. Since classical communication 
also cannot affect entanglement, entangled state classification can likewise involve local operations in con- 
junction with classical communication. A rigorous mathematical characterization of state classification was 
undertaken by Diir, Vidal, and Cirac |42|. Their development was based on states that could be obtained from 
each other through invertible local operators (ILOs), for example, in a three-qubit system, 

|(/)) =A(g)S(g)C|v/-> and |va) = A"^ 05"' (g)C"'|(/)) . (8) 

If this relation holds for two states, they are equivalent under local operations and classical communication 
(LOCC). If this transformation can occur with any finite probability, no matter how small, then the two states 
have equal entanglement under "stochastic" local operations and classical communication (SLOCC). 

For two-qubit states, there exists only one equivalence class of interconvertible entangled pure states, that 
of the Bell states. For three-qubit systems, the SLOCC classification scheme identifies two distinct classes of 
genuinely tripartite entangled pure states. The W class, which retains maximum bipartite entanglement after 
tracing out one of the qubits, is written generically as 

=ai|001>+c72|010)+d4|100) . (9) 

The GHZ class, represented by 

|G/fZS) =c7o|000)+d7|lll) , (10) 

is genuinely tripartite entangled and exhibits no entanglement at the bipartite level. These genuinely tripartite 
entangled classes cannot be converted into one other, even with a finite probability, and are obtainable by 
projections from the generic tripartite state 

|1'3) = ao|000)+a4|100>+a5|101) + a6|110) + a7|lll) (11) 

in J^BC = ■j'lfA®'^^ ®J^c, where the complex coefficients satisfy ^, = 1. 

Although many classifications exist for the multipartite system, studies in ESD and BNSD have primarily 
focused on the generalized «-qubit W class and GHZ class of states. The generalized W-class is represented 
by 

\Wn) = ^(|100...00) + |010... 00) + .. . + 1000...01)) , (12) 

and is characterized by its retaining maximal entanglement under the tracing out of a single qubit as well as 
exhibiting genuine n-partite entanglement. The generalized GHZ-class is represented by 



\GHZn) = ^(|000... 00) + 1111...11)) , 



(13) 
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and is, in contrast, characterized by its genuinely n-partite entanglement, but exhibiting no entanglement when 
any qubit is traced over. 

The classification of pure three-qubit states has been extended to the classification of mixed three-qubit 
states |3|. The generic class in Eq.fTTImay be transformed with a finite probability through local opera- 
tions into the GHZ class of states \GHZ^J^^The GHZ states may themselves be transformed into the other tri- 
partite classes: GHZ, W, biseparable (B), and separable (S) states through stochastic positive-operator-valued 
measures (POVMs) in the following hierarchy 1421 . S C B C W C GHZ. Although there exist some classes 
of multipartite entangled mixed states, such as a generalized n-qubit Werner state composed of a maximally 
mixed state and an n-qubit GHZ-state, these have yet to be considered in the discussion of ESD from noise 
acting on a generalized multipartite state. 

All these entangled-state classifications are based on the requirement that local operations do not increase 
the entanglement of a state because it is global in character. This makes the discovery of ESD and BNSD 
under local noise, although clearly not in conflict with this requirement, somewhat remarkable because a 
global property is seen to be destroyed by local influences alone. Indeed, even very weak local noises are 
capable of destroying these local properties when acting on a range of initially entangled states. 



3.1.2 Entanglement measures 

For studies involving two-qubit entanglement, concurrence C(p) and the closely related entanglement of 
formation Ef{p) are the measures most often used, because they are valid for both pure and mixed states. For 
a two-qubit density matrix p^g, the concurrence is 

— — — (14) 



C(pAB)=max 0,^/^^- \/%- \/?l4 



where the argument of the concurrence function, A = vAi — VA2 — V^4^ is a function of the eigenval- 

ues A, (i = 1,2,3,4), ordered by decreasing magnitude, of the matrix, Pab = Par (c^ ® CT^) Pab (^y^® ^3?)' 

where p^g is the complex conjugate of Pab, o^^^^ the standard Pauli matrix acting on qubit A(B) 148111031 . 
and noting Cab = C'(Pab) to simplify notation for later use; the canonical measure of entanglement, the 
entanglement of formation, can be written in terms of the concurrence as follows. 



EfiP. 



1 -r2 

^ "-ab 



(15) 



where h(x) = —x\og2X— (1 — ji:)log2(l —x). 

The entanglement of both 2x2 and 2x3 complex-dimensional systems is based on the Peres-Horodecki 
criterion: if the partial transposition of the joint system density matrix with respect to a subsystem has one 
or more negative eigenvalues, it is entangled 07 811511 . This qualitative criterion can be quantified via the 
negativity, ,yV{p), given by 

-^(Pab) = I|A;^*"^|, (16) 

k 

where \ are the negative eigenvalues of the partial transpose of Pab with respect to subsystem A, taken to 
be the lower-dimensionality one in the 2 x 3 case 199115011 . Alternatively, but also based on the Peres-Horodecki 
criterion, Huang and Zhu have derived necessary and sufficient conditions for ESD based on a principle minor 
method and the negativity of eigenvalues |54 55 |. 

The existence of entanglement measures that are valid for both pure and mixed states is necessary to pro- 
ceed to larger numbers of dimensions or subsystems. For systems beyond the two-qubit and qubit-qutrit cases, 
such measures are known to exist so far only in special circumstances, for example, when extra symmetries 
are present. For this reason, studies of ESD are not easily extendible to larger Hilbert spaces. Furthermore, 
entanglement must be understood as differing from non-locality. Already in the two-qubit system, there exist 
situations in which entanglement and non-locality differ, the Werner states can exhibit entanglement without 
violating a Bell inequality. The closest larger case for a bipartite system is the symmetric 3x3 -dimensional 
case. Although there is no generalized entanglement measure known to exist so far for arbitrary mixed-state 
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two-qutrit entanglement. Caves and Milburn 1221 derived the separability condition for a two-qutrit Wemer- 
Uke state, p^, composed of the maximally mixed component I9/9, and a maximally entangled component 
|f) = (|ll) + |22) + |33))/y3, 

Ps=^-^-^h + e\W){W\, (17) 

with < e < 1 . This state is separable, that is, not entangled, if and only if e < 1 /4. 

Although a mixed- state entanglement measure exists for the next closest symmetric dimensional incre- 
ment beyond the qubit-qutrit case, a general entanglement measure has not yet so far been shown to exist 
for bipartite states of yet larger dimensions for arbitrary mixed states. However, a measure derived by Terhal 
and Vollbrecht |93 1, which included an ansatz later rigorously proven in [43], does exist for the special case 
of d X d isotropic states, those invariant under transformations of the form U ®U* , where U is unitary II49I . 
The d X (i-dimensional isotropic states piso{d) consist of a convex combination of a maximally mixed density 
{d-^)l^2 and a maximally entangled density P{\W{d))) = \W{d)){W{d) |, that can be written as 

PM = ( ^) 1,2 + ( /"(If (^))) , (18) 

where d>2, 1^2 is the d^ x d^ identity matrix, |f (^/)) = (l/y^)lti fidelity F{pi,^id),Pi\W{d)))) = 
\i{pi^o{d)P{\^id)))) compare with f64l, with < F{pi^oid),Pi\^{d)))) < 1 in the definition of isotropic 
states \49\. piso{d) is separable if and only if F (piso(t/),/'(|f ))) < ^criticai('^) = The entanglement of 
formation for the isotropic states Piso{d) for d >2, 



EfiPiso) = ' 



F < 

Ri.d-i{F), Fe 

^^^^(F-i)+logJ, Fe 



I 

d ' 
1^ 4(rf-l 
d^ 

4{d-l) , 



rf2 



(19) 



where (F) =H2 {F)) + [l-^ (F)] log2(^/- 1), ^2 W = -xlogjCx) - (1 - x) log2(l -x), and ^ (F) = 

2 



F+^{d-l){l-F) 



Other more general quantities relating to entanglement exist for large finite dimensions bipartite states, 
such as the generalized concurrence and entanglement of formation derived in II26II27II . These measures give 
a lower bound for entanglement content, but are often difficult to make use of due to the amount of compu- 
tation involved in performing the minimization and maximization required. For increasing dimensions of the 
discrete-variable dx d system, one approaches bipartite continuous-variable systems of infinite dimensions, 
which are well-described by the Wigner function in terms of momentum pi and position Xt, 

W{px,XuP2,X2)=Y,PiWf'{puXi)W]'{p2,X2) . (20) 

In that limit, the continuous-variable analogue of the Peres-Horodecki condition is that, upon mapping p2 
—p2, the state remains entangled if the Wigner function remains a Wigner function. This criterion was proven 
to be necessary and sufficient for a Gaussian bipartite system [89|. Other necessary and sufficient conditions 
were derived for the Gaussian bipartite cases as well, such as those based on the variances of EPR-like 
operators [41 1. Some progress towards entanglement measures for multipartite continuous states has been 
made, but only for very limited cases that aren't readily applicable for BSD studies. For that reason, the 
discussion here of entanglement in continuous-variable systems is limited to the bipartite case. 

Returning to the discrete case, for a three-qubit system, pure-state entanglement can be measured by the 
three-tangle Tabc for subsystems A, B, and C, 

Tabc = Ca[bc] - C'ab - Cic , (21) 

where the bipartite concurrences involved are calculated in the usual way f3U\. Although this measure is not 
defined for mixed states, it is useful for the identification of classes of genuinely tripartite entangled of states, 
such as the GHZ and W class. One form for a class of entanglement measures for mixed multipartite quantum 
states was developed in 1761 . which is expressed in terms of all reduced density matrices from a state W: 



Cyv(f ) = 21-^/2 /(2^ - 2) {W\ |f )2 - £ tr p2 . (22) 
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The main difficulty with multipartite entanglement measures is the lack of easily parameterized closed-form 
expressions for arbitrary and mixed states. Furthermore, there exist partition-dependent qualitative types of 
entanglement, for example, different sets of fc-partite divisions within a full «-partite system. Hence, one is led 
instead to considering multipartite non-classical correlations in terms of the non-local of correlations rather 
than entanglement, as has been done to extend the sudden death of global state properties beyond BSD to 
BNSD, which is discussed further below. 



3.2 First discoveries 

Addressing the relationship between decoherence and disentanglement, that is, the nature of the loss of entan- 
glement in relation to the loss of state coherence, was an important step toward the discovery of Entanglement 
Sudden Death (BSD) and Bell non-locality sudden death (BNSD). Intuitively, both of these quantum effects 
are expected to behave very similarly. BSD is the extreme case in which coherence persists asymptotically 
whereas the entanglement is entirely eliminated in finite time. In particular, the discovery that they may decay 
at different rates was a hint that BSD can occur. 

An important early step was taken by Yu and Bberly using an exactly solvable toy model consisting of 
a bipartite two-level system coupled to a purely phase reducing reservoir II107II . Their main discovery was 
that entanglement, which they measured using concurrence, decayed at a different rate from the coherence as 
measured by the reduction of off-diagonal density matrix elements. Specifically, the timescale of disentangle- 
ment was always less than or equal to the timescale of decoherence, f^is < ^dec- Other such studies followed to 
clarify the relationship between the rates of decoherence and disentanglement including [ I08ii95i i8i i57j . Fig. 
[T] shows the relationships between subsystems and their environments for a variety of studies. 




Fig. 1 CoiTelated subsystems Sk each in their own noisy environments Ek (k = 1,2, ...,n}. The region enclosed by the dashed 
line represents scenarios typical of two-qubit ESD studies and may represent, for example, two electrons on separate positions 
of a solid state lattice each subject to localized noise influences. Various studies may stipulate further conditions, including: (i) 
each localized subsystem-environment Sk-Ek is noninteracting and nonlocally separated from the others, with only information, 
entanglement, and nonlocal correlations connecting them, (ii) subsystems may share a common environment that affects them 
collectively, (iii) different sorts of noise that are present, such as quantum, classical, dephasing, or amplitude damping, and (iv) 
existence and strength of couplings between subsystems and environments. 



The historically first discovery of BSD was that of Diosi in a bipartite two-level electron system subjected 
to a classical white dephasing environment of magnetic fields (36\. This model includes depolarizing noise 
with decoherence timescale T, the Pauli spin vector a of 2 x 2 Pauli spin matrices, and system density p. 



dp 
dt 



1 

Tt L 



(23) 



which can be interpreted as the statistical average of the wavefunction obeying the Schrodinger equation, 

dY 



dt 



= -- 0)G XI/, 



(24) 



with ft) representing the random magnetic field producing the noise. Using theorems in 152118311 concerning 
disentangling quantum channels, complete disentanglement was shown by Diosi to occur upon satisfaction of 
some conditions proved true using the following arguments. A nonseparable, that is, entangled state is defined 
as one that cannot be written as 

P = 52 PjPsysJ ® Penv J- , (25) 
j 
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where {/?/}, {Psysj }^ {Penvj } are the set of weights, system states, and environmental states, respectively. A 
linear map ^{t) that time-evolves an initial state in time, p(f) = ^(f)p(O), is called entanglement breaking 
if and only if 

■^Po = Y^PjPj, Pj = tr[Pjp] , (26) 

j 

where the set {p/} are densities independent of the original state po, and the probabilities {pj} are obtained 
from {Pj}, the set of positive operator valued measures (POVMs) satisfying Pj > and Y^jPi = I> where I is 
the identity matrix. The solution to this evolution equation can be written as 



p{t) = ^{t)p{0) = ^ l + e-'/^atT ap(0) 



(27) 



where t is the phase decay timescale and ct is the Pauli vector. ESD takes place if and only if 3e~'/^ < 1 . Diosi 
proved this to be true by a constructive proof ||36]| demonstrating the existence of ESD at a disentanglement 
time 

fdis > Tln3 . (28) 

Diosi's discovery of ESD was extended by Dodd and Halliwell ||39ll and Dodd ||38l in a continuous variable 
system prepared in a Gaussian coherent Einstein-Podolsky-Rosen state in position and momentum, subjected 
to thermal environmental noise and negligible dissipation. The Wigner function W{pi,xi,p2,X2) = VK(zi ,Z2) 
evolving according to 

dW pi dW P2 dW d^W d^W 
at m axi m 0x2 op^ op2 

can be solved for the time- dependent solution 

iy,(zi,Z2) = J dh\dh'2g(^i-z\;A)g(z2-z'2;A)Wl)iz[,x'2) , (30) 

where the matrix 

A = Dt(^J Jl^') (31) 
\t/m2t /im^J 

denotes the case of a dissipation-free free particle. The Wigner function was shown to evolve into the separable 
form 

W{px,Xx,p2,X2)=Y,Pi^t'{PUXl)^f{P2,X2) , (32) 

satisfying the continuous-variable Peres-Horodecki condition at a critical disentanglement time 

^ \ 1/2 

nm ^ ' 



t>0.21^-) , (33) 

where D is the dissipation constant and, therefore, exhibiting ESD. Fig. |2] shows the general relationship 
between asymptotic disentanglement and Entanglement Sudden Death. 

A discrete-variable model of spatially separated atoms in a cavity subjected to vacuum noise leading to 
spontaneous emission was next used in the search for additional examples of ESD |109|. Letting Hat, H^m, 
and //int represent the Hamiltonians for the atoms, cavities, and their interaction, respectively, the Hamiltonian 
for the total system is //tot = H^x + //cav + ^^int> with 

Hat =^-(0KO^+^-m , (34) 

/^cav = '"k4'^k + Vk4^k , (35) 
k k 

^^int = I (gla\l + g^,aU^) + £ (fio^bl + halb^) , (36) 
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Fig. 2 A generic plot of entanglement versus timescale, both in arbitrary units, showing the difference between asympototic 
disentanglement (dashed line) and ESD (solid hne), where the entanglement goes abruptly to zero. In both cases, quantum 
coherence decays exponentially. This plot is representative of the qualitatively different ways the entanglement of a system may 
evolve. 



where {gk,fk} are the system-environment coupUng constants, d{d'^) and S(S^) the cavity field lower(raising) 
operators, (7/ (/ = x.y,^) are the Pauli operators, and o± = ^{ax±i(Jy) are the atomic raising and lowering 
operators. In the two-qubit basis 



|1)ab = | + +)ab , |2)ab = I + -)ab ,|3)ab = 
an important class of mixed states with a>0,d= 1 —a, and h 



)ab , |4)ab — I )ab 

: z = 1, is represented by 



Pit)- 



/a{t) 
b{t) z{t) 

zitydt) 

V d{t) 



(37) 



(38) 



Its evolution given the operator-sum representation is 

4 



(39) 



where the operators representing amplitude damping noise, which satisfy the CPTP relations, can be written 
^i = ('nV)^f^?V^.= f^^?Vf ' 2V (40) 



^^3 



'7a0 
,0 1 


(Oa 



Ts 
1 




(Ok 




cob 1 



(41) 



where yA(B) = yA(B) (^) = characterizes the decay for subsystem A(B), described by the rate parameter 

rX^(B) and a)A(B) = y^l — 7a(b)- ^ Markov approximation, the two subsystems decohere at the same rate, so 

that Ia = Ib = J" with analogous relations, 7a(?) = "fsit) = y{t) and C0a(?) = (Oa{t) = (o{t). The concurrence 
is given by 



c(p(0) 



3 



:{o,y(OV(0} , 



(42) 



with f{t) — I- y/a{l - a + 2(o'^ + (o'^a). The satisfaction of the inequahty, 1 - a( 1 - a -|- 2(0^ + (o'^a) < 0, is 
a sufficient condition for concurrence to be zero. Taking the case where a — I, ESD occurs in the timescale 



?dis= pin 



2 + V2 



(43) 
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Since the two-qubit system is important in many contexts, ESD was further developed in this simple, 
but interesting system. The study of larger dimensional constituent subsystems of the bipartite state then 
followed, starting from the next discrete increments in the bipartite 2x3 and the 3x3 cases, proceeding then 
to the arbitrary finite dimensional d xd case. As progress was made for increasing dimensions in the bipartite 
case, research also proceeded involving multipartite states of systems of an increasing numbers of qubits. 

After the preliminary discovery of ESD using the simplifying assumptions of the Markovian approxi- 
mation, further studies relaxed these conditions for a better undertstanding. Non-Markovian systems were 
considered that involve temperature effects, different sorts and combinations of noise, memory effects, as 
well as for different coupling strengths between various subsets of subsystems and environment. Although 
ESD provides an intriguing avenue for mathematical research, for example, its geometrical interpretation, its 
ultimate importance lies with its physics. To this end, its existence was proved in a wide variety of physical 
contexts such as in cavity QED, quantum optics, electrons on a solid state lattice, in superconducting systems 
like SQUIDs, and even in relativistic contexts. After much theoretical development, experimental studies then 
confirmed ESD. These further developments did not proceed in isolation; there is significant overlap among 
them since ESD is such a generic phenomenon. 



3.3 Theoretical developments 

ESD is still a relatively recently discovered phenomenon, about which results have been theoretical ones. 
These results are now surveyed in order of increasing system complexity. 

3.3.1 Two qubits 

The first paper showing ESD in a two-atom system under quantum vacuum noise leading to spontaneous 
emission 1 109] of Yu and Eberly was followed by one examining the effects of "classical" noise 1 1 10|, that is, 
phase damping in a large class of two-qubit mixed states that often arise in physical contexts and includes the 
pure EPR-Bell states and the Werner mixed states. Two types of noise effects were considered, both global 
and multi-local, with the Hamiltonians 

/^global = - l^^Bit){a^ + a^) , (44) 

/^multi-local = - ^M(^AC7f + /^bC?^) , (45) 

where jU is the gyromagnetic ratio and B{t), bp^{t), and b-s,{t) are the classical Markovian, white noise fields 



induced by an external magnetic field. In the standard two-qubit computational of Eq. 37 the "X-states" are 
those that can be written 



Pab(0 = 



[a 

Q b zO 
z* c 
\w* OOd, 



(46) 



with the usual density matrix conditions such as Hermiticity, positive semi-definiteness, and normalization, 
a + b + c + d = 1. When subject to such classical dephasing noise describable in the operator sum repre- 
sentation similar to their earlier study [1091, exponentially decaying off-diagonal terms appear for both the 
global and multi-local cases. Using concurrence as the entanglement measure, the critical time for complete 
disentanglement due to global dephasing noise was found to be 

1 Iwl 

?criticai= ;r7^1n^ , (47) 

where F parameterizes the strength of dephasing and provided b^O and c 7^ 0. For the multi-local dephasing 
case, disentanglement always occurs. Thus, ESD may occur due to classical dephasing noise and so are 
not confined only to either quantum or amplitude damping noise. The authors emphasized that this effect 
is independent of the specific entanglement measure considered. It is also independent of basis, because 
entanglement is so, once a decomposition of the system into subsystems has been chosen. 
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An important later discovery concerned the non-additivity of weak noise influences in ESD. Yu and Eberly 
considered the effects of weak dephasing noise and weak amplitude damping noise, each independently and 
when combined fTTl |, described by Hamiltonians and matrices similar to those used in their earlier studies of 
ESD 1 109 1 10 1 . The concurrence was examined for the X-state density matrix as in Eq.|46j with the following 
form, parameterized by the real quantity A, chosen for simplicity and without loss of generality. 



/I 0^ 
4 A 
OA 4 

Vo 0, 



(48) 



At initial time, C;l(0) =2A/9. For dephasing noise, fl(0) = fl(f) = l/9,ii(0) = t/(f) = 0,z(0) = A/9 ^ z(f) = 
A/9exp[— 12?], yielding a concurrence of 

cf (f) = (2A/9)exp[-r2f] , (49) 

signifying asymptotic decay. Under amplitude damping noise, one finds thatz(O) = A/9 ^ z(f) = A /9exp[— Fif] 
1 /9 -> a{t) = 1 /9exp[-2rif],t/(0) = ^ d{t) = 1 /9ft)f + 8/9ft)i^ with ©i = y'l -exp[-rif] yielding a con- 
currence of 

2 " 



exp[-rif] , (50) 



also signifying asymptotic decay. In the above two expressions, Fi and representing the strength of dephas- 
ing and amplitude damping, respectively. When both dephasing and amplitude damping noises act simulta- 
neously, their combined influence on the state was expected to be additive also. However, the concurrence is 
given by 

Cf •+^""P (0 = 2max jo.Ae'^^' _ ^o)f + 8ft)2| (51) 

in that case, showing that the independent weak dephasing and amplitude damping noises are not additive, 
because the concurrence may go to zero in finite time and so that there is ESD. This study showed that the 
effects of noise on quantum systems, compared to noise acting on classical systems, are not just quantitatively 
different, but qualitatively so. Yu and Eberly obtained further analytic results for these X-states in a follow-up 
analysis 11131 . which also included state-equalizing noise and considered a subset of the X-states, the Werner 
states, given by 

l-f 4F-1 

Pw=^l4 + ^— |f")Cf-| , (52) 

with 1/4 < F < 1, where \W~) is the Bell singlet state. Among their results, it was shown that the Werner 
states subject to amplitude damping noise exhibited ESD only when equal to or less than a new critical fidelity 
Fcrit ~ 0.714, showing that the state was more robust against amplitude damping noise than to dephasing 
noise. ESD was shown to exist for state-equalizing noise as well. Another point here was that initial "local" 
operations may act to protect from ESD, an important point related to decoherence-free subspaces. 

The X-states were also considerd in an analysis by Ikram, Li, and Zubairy 1561 . who considered a two- 
qubit system interacting with a variety of dissipative noise environments (vacuum, thermal, and squeezed 
reservoirs) and showed that ESD occurs in all cases for squeezed reservoirs, as well as for thermal reservoirs 
with a nonzero number of photons. However, for a vacuum reservoir where the mean thermal photon numbers 
are zero, in = n = 0, either ESD or asymptotic entanglement may occur as a function of initial conditions. 
This particular study, similarly to many other studies of ESD, proceeded along the direction involving the 
consideration of different types of noise environments. However, some authors, have studied a bipartite state 
of two-level atoms subject to a controlled laser-induced Stark shift, for example, Abdel-Aty and Moya-Cessa 
121. Using the negativity to measure entanglement, these authors showed that the Stark effect can produce 
either asymptotic entanglement or ESD. As stated in Sec.|2] alternative methods to study ESD may involve 
"intrinsic" decoherence, in contrast to the extrinsic decoherence from the open quantum systems viewpoint. 

Further analysis of ESD in X-states was performed by Ban [ 14 | from the perspective of the relationship of 
bipartite entanglement, as measured by the concurrence, and the phase correlation in bipartite X-states before 
and after the onset of ESD. The phase correlation function of the X-state Pab is given by 

C(0a,0b)=2(|w| + |z|) , (53) 
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where \w\ and |z| are the off-diagonal elements of the X-state density matrix as in Eq. 46 Ban showed that 
before BSD occurs, 

C(pAB(f < fdis)) < C(0A(f < fdis),fe(f < fdis)) • (54) 

At the time of BSD, 

C(pAB(f = fdis)) =C(0A(f = fdis),0B(f = fdi.s))- (55) 

In other studies of the two-qubit systems, Jamroz considered BSD in a two-qubit system due to sponta- 
neous emission and showed that local unitary operations may indeed change the time for BSD onset from 
finite to infinite 161]. Additionally, it was shown that for a large class of states including Werner states and 
pure states, a state may satisfy the two-qubit CHSH inequalities in finite time, even when BSD does not take 
place and the time for complete disentanglement is infinite. Jordan, Shaji, and Sudarshan have discussed the 
effects of "non-" completely positive maps and their expression as local interactions that may increase or de- 
crease the amount of entanglement present [63|. Furthermore, they discussed how noise influences that each 
separately disentangle asymptotically in time, when combined, may in finite time destroy all entanglement. 
Li, Chalaput, and Paraoanu consider a reservoir damping a two-qubit system with various sorts of noise such 
as: transverse, thermal squeezed, and longitudinal thermal |69 |. The novelty of that study was that they were 
also subject to a continuous driving force that may act to decrease the time for BSD onset. Thus, external 
local operations can both protect from BSD and may hasten it. The development of BSD modification through 
local operations was studied by Rau, Ali, and Alber |80|, who showed that multi-local unitary operations 
may control the onset of entanglement sudden death by speeding it up or slowing it down to the point that it is 
never reached. Furthermore, they discovered that when such operations are performed on only one subsystem, 
it may change the rate of disentanglement, but cannot prevent it. Dajka, Mierzejewski, and Luczka consid- 
ered the control of BSD in "non-" Markovian systems ll35ll . Lastly, Maniscalco, Francica, Zaffino, Gullo, and 
Plastina even considered preserving entanglement and so preventing BSD through the quantum Zeno effect, 
which involves periodic strong projective measurements ||73]| . 



3.3.2 Extensions: Higher dimensions or more subsystems 



The initial study of entanglement sudden death (BSD) in two-qubit systems was in the mean time quickly 
extended and complemented by research concerned with verifying and developing a better understanding of 
BSD in finite-dimensional systems by increasing the number of dimensions and subsystems, eventually ap- 
proaching that of the continuous variables. The existence of BSD in the next possible discrete-dimensional 
increment from the two-qubit system was first performed by Ann and Jaeger |10| for a qubit-qutrit system 
using a methodology analogous to that of Yu and Bberly [ 1 09 1 based on the the operator sum decomposition 
to describe multi-local dephasing noise operators acting on two-qubit systems and later for pairs of identical 
subsystems having any finite dimension. This line of research was carried out with an eye toward demon- 
strating that effects that are pervasive in the universe can infinite time destroy entanglement in a broad range 
of systems first prepared in entangled states, not only those of the very simplest case of pairs of two-level 
systems. 

In these studies, dephasing noise was shown to induce a decay in the off-diagonal elements of the density 
matrix of the system of interest in the form of an exponentially decay term and so decoherence only in the 
hmit of infinite time while entirely destroying entanglement on a firute time scale. In this first example, the 
time-evolved qubit-qutrit density matrix was found to be 



pAB{x,t) 



( \ 










loo 



xy(t) \ 




\ 










\ 
1 



\xy{t) 








1 

4 



(56) 



where the diagonal elements < x < ^ were chosen so that reduced subsystems exhibit no coherence and in 
order to satisfy Hermiticity, unitarity, and positive semi-definiteness, and y(t) = exp[— Ff] are the exponen- 
tially decaying decoherence factors with strength F. The entanglement measured by the negativity for this 
state was found to be 

A^[PAB(jc,f)] =max|o,xy(0 - ^1 . (57) 
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Entanglement was here seen to be completely lost when y(f) = 1/8jc at a time f^is = 8jc/F, a finite time for the 
specified range of x and for F > 0. This results also showed that ESD is not restricted to identical subsystems. 

Considering the next higher dimensional increment, ESD in states of a qutrit pair was demonstrated by 
Che9kinska and Wodkiewicz in t24 25] by examining the two-qutrit generalized Werner state 

p. = ^I^®lVe|f^''>(f^''|, (58) 
< £ < 1 with entangled pure state and the pure entangled state 

|f) = 4?(|l^)®|l'^) + |2A)®|2'^) + |3^)®|3'^)) . (59) 
V3 

In the two-qutrit basis {|00), |01), |02), |10), |11), |12), |20), |21), |22)}, the density matrix evolution is de- 
scribed by the operator-sum decomposition of the time-evolution operator with the following Kraus matrices, 




Vl -e-'^i' 0\ /OOv^ 



e 



Ko= Oe-^ \,Ki=\0 0],K2=\00 |, (60) 

0/ \00 

where Ai and A2 are the Einstein coefficients representing the rate of amplitude damping and the Kraus 
matrices satisfy the CPTP relations. The explicit time-dependent density matrix yields the explicit form of 
Caves and Milburn's separability condition ll22ll . 

, , , x 1 

(61) 
and so 



i(f) = I (^2e2^" + 2e-z^2'+2e-z(^i+^2)+e-^i' + e-^2') < ^ 



s{t = 0) = e<^, (62) 

indicating separability; for the maximum value £ = 1, separability was seen to arise in finite time. 

The demonstration of ESD in a bipartite system of a pair of identical subsystems of arbitrarily large 
finite dimensions was then taken by Ann and Jaeger |9|, by considering the case for which a mixed state 
entanglement measure exists for arbitrary d > 2, namely, the d x d isotropic states. The isotropic states are 
those that are invariant under U ®U* transformations and are of the form 

PM = ( ^) hfl + ( Pmd))) , (63) 



defined in Sec. 3.1.2 This state is separable when F (piso(c/),P(|f ))) < Faitica\(d) = c/~ , according to Terhal 
and Vollbrecht's expression for entanglement of formation discussed in Sec. |3.1.2} Specifically, the isotropic 
states must satisfy "both" of the following conditions, 

(i) F(pi,o(^/),P(|f(f = 0))))>Feri,icai(^/), and (64) 

(ii) F{pUd),P{\W{t<c^))))<F,ntic-Ad) (65) 

for some finite time f, for there to be ESD. These conditions were both shown to be satisfied for an isotropic 
state subject to depolarizing noise; an initially entangled state becomes separable in finite time despite the 
persistence of state coherence. This demonstrated ESD indxd isotropic states for arbitrary finite dimensions 
d > 2. A later study also examined entanglement evolution of isotropic states [94|, but by making use of the 
G-concurrence |46 1, which involves the multiplication of "all" of a state's Schmidt coefficients. Therefore, if a 
single Schmidt coefficient vanishes, the G-concurrence vanishes for a c/-level system. That study specifically 
showed that a d x d isotropic state subject to depolarizing noise, parameterized by a strength F, undergoes 
ESD at the k"' level (2<k<d)ina time tk = ln[{d^ - \)/{dk ~d~ l)]/{2dr). 

As the dimensions of the bipartite case are indefinitely incremented, one approaches the case of infinite 
discrete dimensions and the continuous variables case. The work of Diosi f361, Dodd and Halliwell f391. 



and Dodd Ii38i discussed in Sec. 3.2 on continuous-variable states was among the first to demonstrate ESD 
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under local noise. Further work is still required to understand better and clarify the nuances of the relationship 
between the finite and infinite dimensional cases. 

Investigations of BSD in systems of larger numbers of constituents has also very recently been pursued, 
for example, in a study involving BSD in a three-qubit system composed of three atoms A, B, and C each 
interacting only with their own isolated cavities a, b, and c was recently performed ||75l . The total Hamiltonian 
for the system is composed of three Jaynes-Cummings Hamiltonians, 

H = ^ \^z + <^a^iai + g(a ' o'_ + aa'+ ) , (66) 

i=A.B.C 

where (Oq and (O are the atomic and cavity frequencies, respectively. For initial tripartite GHZ and W states, 
the authors were forced, in fact, only to study negativity for bipartitions of the system: A-BC, a-bc, A-bc, 
Ab-c, aB-C, and a-BC, because negativity is not appropriate for mixed tripartite states. At the bipartite level, 
BSD was shown to exist for this set of bipartitions for the tripartite GHZ state, but not for the W state. Since 
these states exhibit a high degree of symmetry, this study strongly suggests, but does not definitively con- 
firm genuine tripartite BSD, because a mixed-state entanglement measure does not exist, despite the authors' 
claims to have done so. 

The entanglement properties of multipartite GHZ-states and W-states were earlier studied using the spe- 



cific form of the multipartite concurrence in Bq. 22 [20 1. The time evolution was studied using the master 
equation 

do ^ 

= £ (I(g)---(g)I)(g)L^(g)(I(g)---(g)I)p , (67) 

k-l N-k 

where p and pf^ are the full system and reduced density matrices, respectively, and Lj^ are the Lindblad opera- 
tors representing multi-local weak Markovian noise, 

p. 

LkP=Y,= W (2«/Pa] - a] atp ~ pa] a,) , (68) 



where J/ represents the strength of the noise and the creation(annihilation) operators as defined in Sec. 3.2 
The concurrence was computed numerically for initial three-qubit GHZ-states and W-states; it was found that 
despite coherence only asymptotically decaying for the zero-temperature case, entanglement vanishes in finite 
time for non-zero temperatures. 

In general, studies suggest that entanglement under environmental influences is more robust for multi- 
partite W-states than for the GHZ-states for a wide variety of situations such as scaling in the number of 
subsystems, temperature, and time. Multipartite BSD due to multi-local noise was also addressed by Aolita, 
Chaves, Cavalcanti, Acm, and Davidovich 1121 . These investigators considered a class of multipartite states in 
which the entanglement calculation was considered in the most balanced partition, and therefore effectively 
reduced to a four-dimensional problem and quantifiable by the negativity. They found that even though the 
BSD time is an increasing function with the number of qubits, the more important timescale for entanglement 
to be arbitrarily small is instead a monotonically decreasing function with the number of qubits. 

The existence of BSD for ranges of initial states in such a wide variety of contexts strongly suggest that 
BSD is a generic phenomenon in all quantum systems in specific classes of states. Given that, in general, 
systems may enter any of a broad range of states, it appears unlikely that entanglement will be found to 
persist in nature, as intuitively foreseen by Schrodinger in his cat paper l,86J . 



3.3.3 Geometry 

In addition to its extension to a greater number of dimensions, the study of BSD has also been developed from 
a geometric perspective. Bntanglement evolution and BSD in a two-qubit system are examined by Cunha |[33l 
and Drumond and Cunha |40|. State evolutions were shown to have different asymptotic states: (i) when the 
asymptotic state is separable, BSD may exist as a feature of the time evolution; (ii) when the asymptotic state 
exists on the border of the entangled and separable states, both sudden death and asymptotic disentanglement 
may exist; (iii) when the asymptotic state is entangled, in which case entanglement sudden birth may be a 
feature of evolution; and (iv) there may exist a multiple asympototic states, in which case some combination 
of all three of the previous may exist. 
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Paz and Roncaglia considered this geometric interpretation of asymptotic state evolution in an explicit 
physical model by examining two resonant quantum harmonic oscillators coupled to a mutual environment 
using the quantum Brownian motion master equation |77|. Taking into account model parameters including 
temperature, ohmicity, and squeezing factors, they showed the existence of and mapped out the boundaries of 
the three phases in this system: sudden death, sudden death and sudden birth cycles, and no sudden death. 

The degree of entanglement was considered to be the distance between a quantum state and the nearest 
separable state was considered by Shimony |87|, who defined the distance as 



where |S) is a normalized separable state and the minimization is taken over all such states. The notion of 
entanglement distance was extended by Yu and Eberly 111211 . who consider the notion of an entanglement dis- 
tance measuring the distance between a given state and the boundary of separable states with entangled states 
and defined by the argument of the concurrence function A. A > 0, A = 0, and A < represent, respectively, 
the cases where the state is entangled, purely separable, and "both" mixed and separable. For example, for 
states that stay only within a separable or an entangled region, and do not become entangled or lose entan- 
glement, respectively, A (?) does not change sign. In contrast, for those states exhibiting entanglement sudden 
death and birth, A(/) may change sign in time. By using this measure to examine the entanglement time evo- 
lution, a rich understanding from a different perspective is obtained. Fig.[3]provides a visual representation of 
asympototic state evolution. 



Fig. 3 A geometric picture of how initially entangled, non-separable states may evolve in time to final asymptotic states. (A) 
Asymptotic disentanglement: a separable state is approached asymptotically in time. (B) Entanglement sudden death (ESD): 
entanglement is completely and abruptly lost in finite time. (C) Coexistence of multiple asymptotic states, entangled and/or sep- 
arable. (D) Entanglement is maintained for the entire state evolution. (E) Entanglement sudden deaths and births: entanglement 
may abruptly reappear after ESD, furthermore, this may occur in cycles. 

Although such a geometric perspective for two-qubit states is illuminating, the results of this analysis pose 
further questions as to the proper geometric description of other entanglement related concepts such as bound 
entanglement, higher-dimensional entanglement, multipartite entanglement, as well as ESD for nonlinear sys- 
tems. Furthermore, these geometric ideas ought to be describable in terms of previously developed geometric 
characterizations of entanglement, such as that of the Minkowskian length associated with the Stokes tensor 



3.3.4 Non-Markovian dynamics 

Many analyses of Entanglement Sudden Death were done assuming the Markovian approximation, where 
calculations are relatively clean and simple analytical results may be obtained, leaving a wide range of less 
pervasive phenomena to be discovered. In particular, the case of non-Markovian noise was rarely considered 




(69) 




m 
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until recently. When the simplifying assumptions of the Markovian approximation are relaxed, many interest- 
ing phenomena arise that aren't possible under that approximation arise, such as backreactive effects, memory 
effects, and entanglement creation. 

The studies so far performed all proceeded in the same following way. First, the Hamiltonians are specified 
for the system //sys and environment //env, along with the interaction Hamiltonian Hi„i that couples them to 
each other. Second, an initial class of states are considered, often with other important initial conditions 
specified. Third, the solution to the time-evolved state is found. Lastly, entanglement dynamics are analyzed 



according to the measures introduced in Sec. 3.1.2 The reader is referred to the appropriate cited publication 



for specific details. When needed, the total Hamiltonian describing the system-environment dynamics is given 



using the relations defined at the beginning of Sec. 3.2 along with the initial conditions. 



Anastopoulos, Shresta, and Hu considered a two-qubit system consisting of atoms coupled to a common 
electromagnetic field and derived the non-Markovian master equation for each of the atoms [6J. Their ap- 
proach is analogous to the study of II109I . but in the non-Markovian regime. The total Hamiltonian of the 
atom-field system is given by //tot = //sys +//int> with 

//sys = hj^cokdldk + hcooia^a^ + a^a^) (70) 

k 

Hint = hY^g, (^l(e-'''-72a^ + e+*-72a«) +d,{e+'^y2at + e-'^'-yia^)) (71) 



where = A/y^ and A is the system-environment coupling strength. The authors considered two classes of 
states, |A) = ^1 1 1) + v^F^lOO) with < < 1, and |S) = x|+) + ( 1 - jc) | -), with |+) = (|01) + 1 10))/y2 

and |— ) = (|01) — |10))/\/2- It was shown that, for class B states, the results are identical to those for the 
Markovian case. However, for class A states, entanglement behavior is dependent on the distance between the 
two atoms. Furthermore, the entanglement exhibited sudden birth and death cycles. Therefore, we see that the 
entanglement evolution is dependent on initial conditions. 

Bellomo, Franco, and Compagno considered a method to characterize dynamics in the case of subsys- 
tems, each influenced different sorts of local noise 1 16 |. This general model was also applied specifically to a 
non-Markovian two-qubit model in a zero-temperature environment, with the Hamiltonian 

//tot = cooo+a^ + £ cokblh + (a+B + a_B"l") , (72) 



where B = Li^g^bk with g^ being the coupling strength. The remaining studies discussed in this section used a 
similar model. Despite the lack of direct interaction between the two qubits, it is shown that their entanglement 
may exhibit cycles of death and birth. 

Dajka, Mierzejewski, and Luczka 1341 considered the exact non-Markovian dynamics of a separated two- 
qubit system, with just one of them subject to local dephasing noise described by the Hamiltonian is given 
by 

//tot = C7i^ + cjf + £ gk{dl+d)af + £ (Okdldk . (73) 

<:=1 k=l 

The reduced subsystems were analyzed for arbitrary subsystem-environment coupling strengths and environ- 
mental fluctuation frequencies. The following different entanglement behaviors were formal: BSD exists for 
the subohmic and ohmic environments, but there is only asymptotic disentanglement for a zero-temperature 
superohmic environment. In a separate study of a non-Markovian two-qubit system, Cao and Zhen studied 
BSD at zero temperature as a function of coupling strength and they found that if the coupling is sufficiently 
strong, BSD occurs, even at zero temperature 1 19 |. 

Liu and Goan recently studied a continuous variable system in a thermal environment subject to non- 
Markovian dynamics 171]. Two types of joint system-environment coupling models were considered, the first 
one having each oscillator coupled to its own environment and the second both oscillator coupled to a common 
environment. This model is analogous to 1 108| for multi-local noise and collective noise, but using oscillators 
instead of qubits. Bntanglement was quantified by the logarithmic of the negativity. Unsurprisingly, behavior 
was found to be dependent on many factors including the initial state, inter-subsystem coupling, subsystem- 
environment coupling, and the influence of multi-local or global noise. BSD was demonstrated to exist for 
some of these situations and not for others. Another analysis of a similar comprehensive nature as this was 
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undertaken by Bellomo, Franco, and Compagno, in which they studied the entanglement dynamics of a two- 
qubit Werner state as a function of purity, initial state entanglement, and the extent of the non-Markovian 
noise character with similar results fTT]. 

Many models exist with even more complex entanglement evolution characteristics, for example, com- 
plete loss of entanglement near quantum phase transitions |47 |. Although this case would not strictly fit the 
classification of BSD, which depends on time, it may conform to the general idea of BSD if the quantum 
phase transition itself could be parameterized as a function of time. 



3.3.5 Energy and temperature 



Several studies of BSD have also been carried out involving its relationship to energy and temperature, pro- 
viding another perspective on the phenomenon. 

Cui, Li, and Yi considered a two-qubit system coupled to an electromagnetic field and analyzed the rela- 
tionship between entanglement evolution and energy transfer when subject to both open systems and closed 
systems dynamics 021 . Considering the (open systems) dynamics involving amplitude damping and dephas- 
ing, they obtained results confirming previous studies on BSD 1 1091 . The Hamiltonian 

(O e 

H=^{ol + o^) + ^-alai, (74) 

was used to model closed systems dynamics. They found close relationships between entanglement and en- 
ergy transfer, even suggesting that entanglement, expressed by the concurrence, could be a function of the 
energy transfer, C[p(f)] = max {0, /(£)}, and that entanglement exists above a critical energy iJcrit but not 
below it. Furthermore, energy and system evolution may be different after the onset of BSD. A similar study 
was performed by Li, Fu, and Liang [70 1 and showed that an BPR state of two-atoms in a cavity disentangles 
faster for a larger initial mean photon number n and that the maximum concurrence corresponds to minimum 
energy of the system. However, entanglement and energy are neither mapped in a one-to-one fashion nor 
evolve at the same rate. 

The quantitative relationship between disentanglement and system-environment energy transfer has been 
developed in more detail by Yu [1061. A two-qubit system in the form of the X-state in Bq. 46 with w = 
evolves according to the Hamiltonian //tot = ^^sys + ^^im + ^^env. with 

H,y, = ]^Ef,a^+^-E^a^ , (75) 
Wenv =Y^O)xa[ax+Y^Vx b\ , and (76) 

Hm = L(/1(7a4 +Aolai) + l^{8laAbl +gxolbx) . (11) 

X X 

An explicit upper bound for the amount of energy transfer required for complete disentanglement of the initial 
state was found: 

£a(2«+1) 

l"^^' - 2[(2n+l)2 + 2n(n + l)] ' ^^^^ 
where n is the mean number of thermal quanta in the reservoir. 

An avenue of research differing from all other open quantum systems studies involving extrinsic decoher- 
ence was undertaken by Silman, Machnes, Shnider, Horwitz, and Belenkiy, where "intrinsic" decoherence is 
considered by the addition of a stochastic term to the Schrodinger equation [88[. The evolution of the initial 
singlet state |f) = (|01) — |10))/v^ was modeled according to the Schrodinger equation in the stochastic 
reduction framework in the energy basis, 

6\x^{t)) = -m\^,(t))At-\ ^ Q{H,-H{t)f\^,{t))At^\ £ C/(^-H(0)IV^(0)dW,(0, (79) 

where |f (0) is the normalized state vector, W{t) is a Wiener process, is a parameter specifying the reduc- 
tion timescale, and dWi{t)dWj{t) = dijdt. After averaging over the noise fields, the state density matrix was 
obtained by the Lindblad equation, 

U[p{t)] = -i[H,E[p{t)]]- £ ^[H^,[H,,E{pm] ■ (80) 

i=A.B ° 
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Using the negativity to quantify entanglement, they determined the conditions for either ESD or asymptotic 
disentanglement to take place as a function of such model parameters as the stochastic reduction rate and the 
direction of the measurement device relative to the motional direction of the particles of the singlet. 

Many studies of the sudden death of global state properties not only assumed the Markovian approxima- 
tion, but also zero temperature. A natural question is how to extend the results of these studies to the case of 
finite temperature T > 0. In order to address this quesetion, Jakobczyk and Jamroz considered a two-qubit 
system coupled to a finite-temperature bath evolving according to the following master equation, 

^ = lm[o+,pa^] + [a+p,a-]}+^-m[a^,pa+] + [a^p,a+]} , (81) 

where 1] = 7o«(ft)o), I] = 70 (1 +n(ft)o)), n(coo) = ^ , and P = j II6OI . All initially entangled states in a 
noisy cavity were found to disentangle on a finite time scale that is an increasing function of the initial state 
entanglement. Furthermore, they showed the onset of loss of standard Bell inequality violation occurs on a 
shorter timescale. Al-Qasimi and James in a similar analysis considered a two-qubit system initially prepared 
in the X-states and coupled to some finite temperature exhibit ESD |5 1. Notably, these classes exhibiting ESD 
contain states that were robust against ESD in the zero temperature regime. 

These studies suggest that ESD is not only a novel and counterintuitive result, but it is more likely than 
long-lived entanglement in typical contexts, which involve energy transfer and temperature effects. Practical 
tasks such as quantum error correction protocols necessarily need also to account for energy transfer and 
temperature-induced effects that degrade entanglement and coherence themselves in addition to the other 
sorts of phase and amplitude damping errors. 

3.3.6 Baths: spins and oscillators 

Many studies of ESD consider a central system of interest composed of spins and/or oscillators surrounded 
by and interacting with a bath of spins and/or oscillators. In this section, the case of spins is mainly con- 
sidered. The dynamics for the central system are then solved for, either before or after tracing over the en- 
vironment variables of the state, at which point the coherence and entanglement dynamics are examined. 
They also usually involve changing subsystem-environment coupling strengths or distance scales and other 
aspects of the model environment. There is often reference to different environments according to their spec- 
tral properties. It is therefore useful to recall that the environmental spectral function in the thermodynamic 
limit is /(ft)) = Aft)''*''^ exp[— ft)/ft)|;rit], where (Ocdt is the environmental cutoff frequency and A is the qubit- 
environment coupling strength; the sub-ohmic, ohmic, and super-ohmic case corresponds to jj. G (—1,0), 
;U = 0, and ji G (0,°°), respectively {cf. [68|). 

Roszak and Machnikowski study ESD in a two-qubit system subject to pure dephasing influences from 
a super-ohmic environment and the effect of inter-qubit distance lISTIl . The two-qubits and environment 

represent excitonic qubits and phonons, respectively. For the initial states Iv^g^') = 1^') ^^d 

IVo^') = ^^^^^72^' entanglement dynamics were analyzed using concurrence. ESD was shown to exist 
as a function of different initial states and temperatures. Additionally, inter-qubit distance was shown to be 
important since ESD only occurs for spatially separated states. 

The importance of distance, mentioned by Schrodinger |86|, between subsystems in relation to entangle- 
ment was carefully considered by Cormick and Paz, who studied a two-qubit system locally coupled to an 
environmental XY-spin chain [31 1. After obtaining an exact solution of the system of interest, they analyze 
how decoherence and disentanglement behave as a function of distance and coupling strength. For strong 
system-environment couplings, quasiperiodic entanglement sudden deaths and sudden births as a function 
of the inter-qubit distance were found. Jing, Lu, and Yang [62| consider a two-qubit Bell state where one 
qubit is not coupled to the environment and the other interacts with a Heisenberg XY spin bath in thermal 
equilibrium. In one case, the bath contains a finite number of spins (A' = 40) and in another case an infinite 
number of spins (i.e. A'^ 00). After finding the time-evolved state using numerical methods based on II37I 
[53 1 and obtaining the reduced density matrices for the two-qubit system by tracing over the environmental 
variables, Psys(0 = trenv[p(0]' they calculated the concurrences for the initial Bell states and found that (i) 
the entanglement evolution was independent of the initial state and whether the spin-bath was composed of 
a finite or infinite number of spins, (ii) a smaller anistropy parameter can more easily allow the system to 
become entangled, (iii) finite temperatures allows for ESD, and (iv) large intra-subsystem coupling can help 
maintain a high degree of entanglement. 
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Ma, Wang, and Cao considered three-qubit systems in a quantum-critical environment consisting of an 
Ising model in a transverse field |72|. They examined the evolution of the entanglement, quantified by the 
negativity, as a function of system-environment coupling strength, environmental degrees of freedom, trans- 
verse field strength, and in relation to state symmetry. The three-qubit states considered include: the standard 
GHZ state \GHZ) = ^(|000) + |111)), the W state \W) = ^(|001) + |010) + |100)), and the three-qubit 

Werner state Ps-Wemer = fis + (1 — p)|GHZ)(GHZ|. For these, they considered the entanglement between 
for the following bipartitions using the negativity: A^ab-c. ^ac-b. ^bc-A; ^a-b. ^a-c. ^b-c> and confirmed 
the existence of BSD relative to these bipartitions as they explored the conditions upon which BSD takes place 
and found that disentanglement can be enhanced by a quantum phase transition when the system is coupled 
weakly to the environment. Although the lack of entanglement in all bipartitions strongly suggests the lack 
of entanglement at the tripartite level, this is not necessarily so; demonstrating this requires a rigorously de- 
fined genuinely tripartite mixed state entanglement measure, which has yet to be found. Sun, Wang, and Sun 
considered a similar model II91I . For the two-qubit system with entanglement evolution quantified by the con- 
currence, the initial pure state | ) = ( 1 00) + 1 1 1 ) ) asymptotically disentangles, whereas the initial mixed 

two-qubit Werner state P2-Wemer ~ P 

|<J>+)(<J>+| + i^l4 exhibits BSD. Similarly, the two-qutrit system with 
entanglement evolution quantified by the negativity, the initial pure state | f ) = ^ ( 1 00) + 1 1 1 ) + 1 22) ) exhibits 
asympototic disentanglement, either to zero or to some constant, and the initial mixed three-qubit Werner state 
P3-Werner = p\^) (f| + ^^^9 bccomes Completely disentangled in finite time. For both the two-qubit and two- 
qutrit systems, the specific rate of asymptotic disentanglement and the actual time of BSD in this situation is 
determined by several factors, including the transverse field strength and the system-environmental coupling 
factors. Nonetheless, the authors demonstrated the basic qualitative result that BSD occurs in this model, 
confirming previous results. 

Lai, Hung, Mou, and Chen consider two-qubits coupled to an XXZ spin chain II66I . The novel feature of 
their analysis was that it goes beyond the usual method involving a central system coupled to the environment, 
and instead uses time-dependent density renormalization group theory (t-DMRG) to non-perturbatively model 
spin-bath dynamics in order to determine the decoherence and disentanglement properties. It was shown in 
this way that the phase of the environmental spin bath is important to entanglement dynamics. However, 
only in the anti-ferromagnetic and paramagentic phases was BSD shown to exist. Using the concurrence, 
the authors determined the time for BSD for spin-bath phase diagram. Numerous other results were also 
discovered, demonstrating the relevance of quantum phase transitions and BSD. 

An important example of the analogous case to spins, that of harmonic oscillator systems, is the a com- 
prehensive analysis of Chou, Yu, and Hu 1281 . who demonstrated that BSD exists under a wide variety of 
situations. 



3.3.7 Other realizations ofESD 

The study of BSD in abstract models was further developed in more physically realistic contexts with an 
eye towards experimental observation and practical application. An example from each of these contexts is 
provided in this section. 

Cavity quantum electrodynamics (QBD), which has been well-developed theoretically and experimentally, 
is well-suited to investigate BSD. One particular cavity QBD model is frequently used in BSD studies - that 
involving two atoms in independent cavities evolving in time under a double Jaynes-Cummings Hamiltonian, 

Hjc = yai^ + + a^a) + (oa^'a+ ^af + g{b^'a^ + a^b) + (ob^b , (82) 

with h = I, (0 the cavity frequency, (Oq the atomic frequency, and g the atom-cavity coupling strength. This 
basic Hamiltonian can be modified in a variety of ways, such as changing the atom-cavity coupling strengths, 
specifying various combinations of atomic and cavity frequencies, and considering different initial entangled 
states. Using this Hamiltonian, Yonaf, Yu, and Bberly |105| considered a four-qubit model with a cavity 
QBD interpretation, such that all reduced density matrices are of the X-state form, as in Bq. |46] with z — 
0. Six concurrences were calculated (Cab. Qb; Cah, C^b, CAh, and Cfia), allowing them to study BSD and 
entanglement revivals, with the stipulation that there is no contact between either of the two atoms or the 
two cavities. Bntanglement evolution was considered a pure exchange and information-based phenomenon 
independent of actual physical interaction. 
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Quantum optics also provides many opportunities in which to better understand ESD. The study of co- 
herence and entanglement in this field has been especially well researched, due to their relation to practical 
quantum communications technologies. For example. Gong, Zhang, Dong, Niu, Huang, and Guo considered 
states of single- and double-photon polarization states subject to local phase noise 1451 . in the form of Eq. 
[46) again with z = 0. The decoherence and disentanglement properties were examined as a function of differ- 
ent frequency spectrum envelopes. They discovered that some colored frequency spectra lead to asymptotic 
decoherence and disentanglement and that others can lead to periodic sudden deaths and births of coherence 
and entanglement due to non-Markovian effects. 

In addition to the discovery and analysis of ESD in cavity QED and quantum optical setups, ESD has also 
been determined to exist in other physical systems. For example, Abdel-Aty showed the existence of ESD in 
a Josephson junction by modeling a two-qubit maximally entangled mixed state, with each qubit of a Cooper 
pair box connected to a reservoir through a Josephson junction 1 1]. Roszak, Machnikowski, and Jacak showed 
that ESD exists on in a system composed of excitons on quantum dots subject to phonon dephasing noise and 
depends on the distance between these two excitons due to the finite phonon phase velocity II82I . 



3.4 ESD: Empirical evidence and experimental proposals 

Given that entanglement sudden death has been shown to exist in the theoretical context, both in abstract treat- 
ments and in realistic models of specific physical systems, it is not surprising that specific experimental tests 
have been proposed and that some of the above results have already been experimentally confirmed. The ex- 
perimental observation of entanglement sudden death in physical systems necessarily requires measurement, 
often requiring resource-intensive full tomographic state reconstruction. This is so because the construction 
of entanglement measures may involve operations that do not necessarily have one-to-one correspondence to 
the physical world, for example, the cutoff function used to compute the concurrence and partial transpose 
operation used in the negativity are abstract mathematical operations. 

Theoretical support for the measurement, characterization, and dynamically monitoring of entanglement 
evolution has been provided, for example, Santos, Milman, Davidovich, and Zagury proposed the use of 
entanglement witnesses to signal the existence or absence of entanglement by monitoring a single measur- 
able quantity that remains constant throughout the evolution of the quantum system 1 85] . Furthermore, they 
proposed an experimental setup that may be used for cavity quantum electrodynamics and for trapped ions. 
CarvaUio, Busse, Brodier, Viviescas, and Buchleitner [21 J generalized the idea of entanglement measurement 
and showed that for experimentally relevant mixed quantum systems evolving under open quantum systems, 
entanglement monitoring may proceed in an "optimal" way based on quantum jump operators. 

Experimental evidence for ESD in a variety of physical contexts such optical setups and atomic ensem- 
bles has already been given. Almeida, de Melo, Hor-Myell, Salles, Walborn, Souto Ribeiro, and L. Davi- 
dovich experimentally confirmed the existence of entanglement sudden death for a two-qubit system due to 
multi-local dephasing and amplitude damping noise in an optical experimental setup involving a Sagnac-like 
interferometer [4]. In the two-qubit system, one qubit is denoted by the horizontal and vertical polariza- 
tions of a photon, the other qubit is the ground and excited state of an atom, and the environment acting 
upon these two-qubit system is the momentum of the photon. The general photon polarization Bell-state 
\'P) = \ a\\HH) + |j3| exp(i5) \VV), where H is horizontal polarization and V is the vertical, is considered and 



may be viewed as a subclass of the X-states, with w = z = in Eq. 46 They consider two initial states: 
defined by |j3p = |o:|^/3 and defined by |j3p = 3|ap. These initial states both have a concurrence of 
C(p) ~0.8 and similar purity, respectively, w 0.91 and ^// w 0.97. |a|, |j3|, and 5 are modified physically 
by a combination of quarter- and half-wave plates put in the pump beam path. The evolution of the system for 
an amplitude damping channel is given by the following map, 

\H)<E)\b) ^\H)iS)\b) (83) 
\V) ® \a) /r^|y) ® \a) + ^\H) ® 1^) , (84) 

where denote orthogonal spatial modes. Entanglement is measured by the concurrence, given by 

C(p)=max {0,2(1 -p)|j3|(|a|-p|/3|)} . (85) 

Here, for the case of |j3 1 < |a|, there is no entanglement if p= 1. In contrast, for |j3 1 > |a|, finite-time disentan- 
glement occurs for /) = |a/j3|. They showed that \ ^fi) undergoes asymptotic disentanglement, with complete 
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disentanglement occuring only when p = I, the case where each individual subsystem is completely deco- 
hered. For Ixjfii), the concurrence goes to zero for p < 1, showing BSD. The dephasing channel is described 
by the following map, 

\H) ® \b) -> \H) (g) \b) (86) 
\V) ® \a) Vl-p\V) ® \a) + ^\V) (g) \b) . (87) 

Here, both states | \ffj) and !<////) exhibit the same behavior; only for p=l do they completely disentangle. This 
study was the first experimental confirmation of BSD. Salles, de Melo, Almeida, Hor-Meyll, Walborn, Souto 
Ribeiro, and Davidovic later expanded on this discovery by a comprehensive analysis of quantum optical 
experiments exploring further BSD for the amplitude damping channel [841. 

Laurat, Choi, Deng, Chou, and Kimble studied the heralded entanglement between two qubits, each repre- 
sented by collective excitations of a cloud of cesium atoms II67I . For purposes of determining the entanglement 
content, the states of the two atomic ensembles are each locally mapped to photonic states. Since entangle- 
ment is a non-increasing quantity under local operations, this mapping process ensures that the entanglement 
in the photonic states cannot be more than that in the atomic ensembles. The relevant density matrix is of the 
X-state form of Bq. 46 with w = 0, given in the photon number ({«,/«} = {0, 1 }) basis, \n) \in), 



/POQ 
poi z 
z* pw 

V pn 



(88) 



with P = pqq + poi + pio + Pn the normalization factor, pij is the probability of / and j ({n,m} = {0, 1}) 
photons existing in each of the photonic modes. The concurrence is then given by 

C(p)=max[0,Co] with Co='^{2\z\^2^mm) . (89) 

Concurrence as a function of storage time begins above zero and goes abruptly to zero in finite time, proving 
the existence of BSD in a system composed of entangled atoms ensembles. 



4 Non-locality Sudden Death 

A significant obstacle for the investigation of Bntanglement Sudden Death (BSD) for mixed multipartite sys- 
tems has been the lack of a generalized entanglement measure. However, the investigation of the sudden death 
of the somewhat similar global state property of Bell non-locality has been feasible. In particular, the finite- 
time loss of Bell-non-locality under noise-induced asymptotic state decoherence. Bell non-locality Sudden 
Death (BNSD), has been theoretically demonstrated, extending research on the destruction of global state 
properties by local noise to three-qubit and larger systems. This avenue of research may also help illuminate 
the relationship between entanglement and non-locality, which have differences the origin of which are not 
well-understood, even for a two-qubit system. For example, Werner showed that a class of mixed bipartite 
two-qubit states can be entangled despite not violating a Bell inequality 1 101 1, that is, some entangled states 
may admit a local hidden variables model [96 1. Similarly, it is conceivable that the study of BSD and BNSD 
could provide a better understanding of other fundamental questions, such as the reason for the specific value 
of the Tsirelson bound ll97ll98l . 



4.1 Non-locality measures 

Bell Non-locality Sudden Death (BNSD) has so far been studied in two- and three-qubit systems. For the bi- 
partite two-level system, there exists a single necessary and sufficient condition for Bell non-locality, allowing 
it to be readily studied. Specifically, non-locality is detected by violation of the Clauser-Horne-Shimony-Holt 
inequality [29 J : 

^i=\ [MaMb +MaM^ +M'p^Mb - M'p^M'^] , (90) 
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where the operators M and M' act in different directions, \^2\ < 2 indicates the existence of a local classical 
model and 2 < \^2\ < 2^/2 denotes the presence of quantum correlations. The upper bound of 2\/2, the 
Tsirelson bound, specifies the maximum reached by quantum mechanics; the algebraic maximum at \^2 \ = 4 
denotes maximal non-locality, with the region 2^/2 < |^2| < 4 containing non-physical "super-quantum" 
correlations that are implied when only a subset of known physical principles are enforced. 

For a three-qubit system, both the classification of initially entangled states and the set of Bell inequalities 
used to detect their non-locality become more complex. (The classification methodology and tripartite SLOCC 



classes are introduced in Sec. 3.1.1 ) For such systems, tripartite Bell inequalities fall into two categories 
||23 |. The first category, the Werner- Wolf -Zukowski-Brukner (WWZB) inequalities 1102111141 . distinguishes 
between tripartite states that are describable by a local classical model and those that exhibit "any" non- 
locality at all. The second category, the Svetlichny inequalities, distinguish states that can be described by a 
hybrid local-non-local model and those that cannot. Let us discuss these in turn. 

For tripartite states to be describable by a local classical model, the "entire" set of Werner- Wolf -Zukowski- 
Brukner inequalities must be satisfied I1102II1141 . For the tripartite case, there exist five classes of inequalities 
encompassing thirty-two total inequalities, with the inequalities contained within each class equivalent upon 
symmetries. Because the conjunctions of these inequalities are a necessary and sufficient condition for non- 
locality, the satisfaction of "all" these inequalities is necessary and sufficient to demonstrate a local classical 
model, however, only a violation of "one" is sufficient to show the existence of non-locality. They are: 

(PI) ^Pi = 2MaMbMc , 

(P2) ^P2 = ^ ( -MaMbMc + MaMbM'c + Mp^M'^Mc + Mp^M'^M'c 

+ M^MbMc +M'pMYiM'c +M'pM'^Mc +M'fJ4'^M'c) , 
(P3) ^P3 = [Ma (Mb +M^) +M'^iMB - M{,)]Mc , 
(P4) ^P4 = MaMb (Mc + Mc ) - MaMb (Mc -M'c) , 

(P5 ) ^P5 = MaMbM^ + MaM^jMc + M^MbMc - M'^^M^M'^ , (91) 

where M, (;' =A,B,C) are projective measurement operators acting on the qubit and the prime denotes an 
alternate measurement. All of these must be satisfied for a state to be describable by a local classical model, 
I (^p/)p I < 2 (/ = 1, 2, 3, 4, 5). This set includes the often considered sub-class of tripartite Bell inequalities 
174111 3 ifol . namely those involving ^P5. 

Although the WWZB set of Bell inequalities provide necessary and sufficient conditions for determining 
whether a state can be described by a local classical model, they cannot distinguish between whether a state 
is describable by a truly non-local model or a hybrid local-non-local model. For this, one needs to consider 
the class of Svetlichny inequalities, represented by 

S = MaMbMc +MaMbMc +MaMbMc +MaMbMc 
- MaMbMc - MaMbMc - MaMbMc - MaMbMc , (92) 

where the measurement operators M(M') are defined as before, |S| > 4 denotes genuine tripartite Bell non- 
locality, |S| = 4^/2 is the maximum attainable through quantum mechanical correlations, and |S| = 8 is the 
algebraic maximum |92|. (Note that there exists another Svetlichny operator expectation 5' that may contain 
have a different phase, but it is equivalent to this one when only the magnitude of the quantity is considered.) 



4.2 Specific BNSD studies 

Ann and Jaeger |58 1 were first to consider the finite-time loss of Bell-inequality violation (BNSD) in a three- 
qubit system. They considered a pure generalized W state, 

=ai|001)+a2|010> + c74|100) , (93) 

with coefficients a, G C (;'= 1,2,4) satisfying Y.i'^f = 1- Using the operator-sum decomposition, the time 
evolution was considered in an analogous way to previous BSD studies, initial state is considered subject 
to multi-local dephasing noise, which affects the state density matrix by the appearance of multiplicative 
exponential decay factors y{t) = on the off-diagonal elements, where F denotes the strength of local 
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dephasing. For an asymptotically decohering state, the finite-time loss of Bell-non-locality is demonstrated 
when "both" of the following conditions involving the two different (non-)locality regimes are satisfied: 

(i) (■^2)p(o)>2 (94) 

(ii) (^2)p(,)<2, (95) 

with f2 < f < °°, f2 denoting the time a violation ceases to exist. The specific measurement operators used for 
qubits A, B, and C were 

Ma = a,®I(g)I , (96) 

m;^ = a,(g)I(g)I , (97) 

Mb = I(8)[cos(0B)c7,-sin(0B)o:v]®I , (98) 

Mb = I(8)[sin(0B)aj + cos(0B)c!:x]®I , (99) 

Mc = I(8)I(8)[cos(0c)cyz-sin(0c)crt] , (100) 

M'c = l(S)l(S)[&in{dc)a^ + co&{dc)ax] , (101) 

where and are Pauh operators, with the rotational angles given by the specific values 9b — ^ and 0c = f- 
For the maximally entangled case where the coefficients of the \Wg) are equal and real, fli = 02 = 04 = 1 / a/3, 
it was shown that complete loss of Bell inequality violation occurs at the finite time scale 

ln(2) 

f2 = ^ . (102) 

In a follow-up publication |fTT|. the authors of jSH improved upon the above initial study of the finite-time 
loss of Bell- inequality violation in W-states by considering the effects of multi-local dephasing noise acting 
on the initial generic class of tripartite entangled states, 

1^3) = do|000)+a4|100)+c75|101)+a6|110)+a7|lll) , (103) 

with a more nuanced examination of non-locality properties using the entire set of WWZB and Svetlichny 
inequalities. Genuinely tripartite non-locality was considered first, using the Svetlichny inequality. Here, the 
maximum quantum violation | (S) | = 4\/2, is reached by the GHZ state with |flo| = |«7 1 = 1/ and all other 
coefficients zero. Under multi-local dephasing noise and for the set of measurement operators in Eq. |101| 
chosen for maximum initial time violation of S, the time-dependent Svetlichny expectation of the three-qubit 
state is given by 

|(5)p(,)|=8\/2|ao||a7k-'^S (104) 

where F parameterizes the strength of the dephasing noise. Tripartite non-locality for the maximally entangled 
GHZ state with do = dj = 1/ a/2 and all other coefficients zero, as distinguished by the Svetlichny inequality, 
is completely lost at a time 

= (.05, 

and for all times afterwards. Figs. |4] and |5] provide a visualization this time evolution. 

Next, the loss of "all" non-locality was considered using the set of WWZB inequalities. The maximum 
quantum violation possible for the set of WWZB inequalities is | (^) | = 4. Determining the non-locality 
properties of using the WWZB inequalities requires examining each of the five independent classes of in- 
equalities, however, the time required for complete loss of Bell-inequality violation was found to be bounded 
from above by the time it takes to be completely lost. Therefore, the relevant quantity that captures Bell 
non-locahty decay of the W state is given by 

|(Sp5)p(„|=8k7o|k77|e-^^^ (106) 

The "entire" set of WWZB inequalities was shown to be completely lost in the finite time scale 

(.07) 
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Fig. 4 At initial time, the Svetlichny expectation (S)p(Q) is plotted as a function of the rotation angles 9b and 9c for the maxi- 
mally entangled tripartite state, \GHZt,) = (|000) + |lll))/\/2. | (S)p(o) I > 4 denotes violation of the Svetlichny inequality and 
genuinely tripartite Bell non-locality. 




Fig. 5 These plots show subsequent time evolution of the Svetlichny expectation as the state is subject to multi-local dephasing 
noise. On the left, | (S)p(,.) | < 4 for all measurement angles 9b and 9c, indicating the lack of genuinely tripartite Bell nonlocality. 
Because this occurs in finite time t^, despite violation at / = and asymptotic decoherence, genuine tripartite Bell non-locality 
sudden death is shown to exist. Subsequent time evolution is shown on the right, for example at f = Sfj > t^, \ (S)p(,) | < 4. The 
set of WWZB inequalities exhibit similar behavior, denoting the loss of "all" generic Bell non-locality in finite time. 

and for all times afterwards. Thus, the tripartite state was shown to be described by a local classical model. 
This analysis used the Svetlichny inequality to demonstrate the loss of genuinely tripartite non-locality and 
used the entire set of WWZB inequalities to show the transition from a non-local state to a local classical 
state; despite state decoherence occuring asymptotically in time, Bell non-locality was completely lost in 
finite time. 

Tripartite BNSD has very recently been applied in other contexts, making use of a similar methodology 
as the initial discovery. Yang, Yang, and Cao considered a system composed of two identical atoms in a lossy 
cavity according to the Tavis-Cummings model II104I . The initial entangled state for the tripartite system is 
composed of the two atoms and the cavity, was taken to be the W state, | V''(0)) = a\^gG) +^|^^'0) +c\gg\), 
where \g) and \e) denote the ground and excited states of the atoms, respectively, |0)(|1)) denotes the two 
different states of the cavity, and + Z?^ + = 1 . The non-locality properties are exmained as a function 
of time for the maximally entangled W state with a = b = c = 1/ They calculated the critical times for 
loss of Bell inequality violation for different relative values of the photon leakage rate K and the coupling 
constant k; t{K = k) > ln(2) / K and t{K = k) = ln(2) /k; when K <k, the tripartite non-locality oscillates 
before decays asympototically. Note that there can exist non-locality revivals at the boundary. 

Qiu, Wang, Su, and Ma similarly considered a tripartite state of spins coupled to an antiferromagnetic 
environment with an applied magnetic field |79|. The Hamiltonian is derived using the spin-wave approxima- 
tion. For an initial separable system composed of the tensor product of an environment in thermal equilbrium 
with a pure tripartite GHZ state, these authors found the time until BNSD as a function the applied magnetic 
field. When the initial system is a W state, there is no BNSD since the state is a decoherence free state in 
that environment. For this reason, they have suggested that the generalized W state is a good candidate for 
quantum information processing purposes. 

Although no experimental studies have yet been performed for the specific purpose of testing the existence 
BNSD, since it was first theoretically explored only in 2008, it is likely that some past experiments designed 
to verify the existence of non-locality itself but failed to definitively do so were experiments in which BNSD 
took place during between state preparation and state measurement. 
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5 Implications and summary 

The above review has considered a wide range of theoretical situations and a number of experimental sit- 
uations in which entanglement sudden death (ESD) and Bell non-locality sudden death (BNSD) have been 
demonstrated. These phenomena have been explored only for the last five years in the first instance and one 
year in the second, leaving considerable room for further study. The existing results already suggest that 
previous studies involving asympototic "decoherence" are readily extendable to the consideration of disen- 
tanglement and the loss of non-locality in finite time in a way that may bear on the foundations in ways that 
the study of decoherence has so far failed to, despite initially high expectations for that earlier related program 
of study. 
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